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ABSTRACT 

The soil water loss during the drying period of any irrigation cycle leads to progressive leaf 
dehydration and consequent rise in canopy temperature. Irrigation scheduling through crop tem- 
perature measurements is a technique well suited for infrared radiometric remote sensing. A soil- 
plant-atmosphere model together with observed hourly weather uata are used to simulate soybean 
(Glycine max L.) leaf water potential, stomata! resistance and canopy temperature at various 
soil water potentials. Using an empirical relation between the stomatal resistance and leaf 
water potential, the model solves for the leaf water potential which satisfies Monteith’s and 
van den Honert’s equations for transpiration. The canopy temperature is then obtained from 
the energy balance equation. The hourly weather data are for several clear sky days during 
summer at Phoenix and Baltimore, and covered a wide range of variables. For soil water 
potentials near field capacity the simulation gives a base line relationship between the canopy- 
air temperature difference (5T) and air vapor pressure deficit (VPD) which agrees well with a 
linear regression equation developed from observations. At high soil water potentials the effect 
of weather variables on the canopy temperature can largely be accounted for in terms of air 
and dew point temperatures. The 5T values at lower soil water potentials are found to be 
uniquely related to the base line 5T values. The simulation indicates that the root zone soil 
water potential may be inferred from observed 5T and the base line 5T calculated from air and 
dew point temperatures. At high soil water potentials the simulated stomatal resistances show 
scatter, the range of which agrees well with a set of observations. The simulated increase of 
stomatal resistance with decreasing soil water potential is also in agreement with a set of ob- 
servations. The simulation somewhat overestimates the leaf water potentials which is suggested 
to be due to some characteristic differences among soybean varieties. 
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SIMULATING SOYBEAN CANOPY TEMPERATURE AS AFFECTED 
BY WEATHER VARIABLES AND SOIL WATER 
POTENTIAL 


INTRODUCTION 

Effective use of irrigation water is rapidly becoming a subject of considerable interest, and a 
method for irrigation scheduling based upon crop temperature measurements has been suggested 
(Aston and van Bavel, 1972; Jackson, 1982). This technique for irrigation scheduling is well suited 
for infrared radiometric remote sensing. It has long been recognized (Tanner, 1963) that crop tem- 
perature can be used as an indicator for soil water stress; as plant available soil water decreases the 
canopy temperature increases. To identify crop water stress more-or-less unambiguously by crop 
temperature measurements, it is necessary to quantify the effects of weather variables on the 
canopy temperature. 

Carlson et al. (1972) observed that air temperature and vapor pressure deficit (VPD) are the 
significant weather variables affecting soybean canopy temperature. For a small range of VPD, 

Ehrler (1973) found a significant linear relationship between the canopy-air temperature difference 
(5T) and VPD for cotton. Using data from different geographic locations, Idso et al. (198 1) showed 
that during a significant portion of daylight periods unique linear relationships exist between 5T and 
VPD for non-stressed alfalfa, soybean and squash. 

The objective of this paper is to develop a soil-plant-atmosphere model for soybeans and study 
the sensitivity of canopy temperature to weather variables, soil water potential, and crop rooting 
density. The model is tested against the observations of Idso et al. (1981) for canopy temperatures, 

Brady et al. (1974) for leaf water potentials, and Brady et al. (1975) and Sivakumar and Shaw 
(1978) for stomatal resistances. 
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MODEL DESCRIPTION 

The Monteith’s (1965) equation for transpiration contains a variable resistance term for 
representing vegetative control of vapor loss. The basic idea of the soil-plant-atmosphere models 
(Federer, 1979; Soer, 1980 and Luxmoore et. al., 1981) is to account for this vegetative con- 
trol of vapor loss by setting up an equation of continuity for the plant water status. If in the 
transpiration process the canopy neither gains nor loses any water then the continuity equation 
can be written as 
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AR n fC n Pa ( e a* “ e a ) / r a . j s ~ jj (0 

L v fA + 7 (1 +_r £< ) 1 * R s + R p 

r a 

where the left hand side is the Monteith’s equation for transpiration and the right hand side 
is the van den Honert’s 0948) equation for water uptake by plant roots. 

The left hand side contains the atmospheric parameters which determine the atmospheric 
evaporative demand and the stomata! control of transpiration. The saturated vapor pressures 
(kPa) at air and dew point temperatures are, respectively, e a * and e a , and their difference is 
the vapor pressure deficit (VPD). The slope of saturated vapor pressure evaluated at approxi- 
mately the air temperature is A (see Jackson et al., 1981). The density (kg m “ 3 and the heat capacity 
(J kg " 1 K" 1 ) of air, are, respectively. p a and C p , 7 is the psychrometric constant (kpa KT 1 ) and L y is 
the latent heat of vaporization per unit volume (J m" 3 ). The aerodynamic resistance r a (s/m) is 
calculated from wind speed U (m/s) at height Z according to Thom and Oliver (1977) as 

4.72 [ In ( ^~ d ) j 2 

r = 2 ( 2 ) 

a 1 + 0.54 U 1 

where the roughness height Z Q and the zero-plane displacement d pan be calculated from the 
crop height CH as (see Bailey and Davies. 1981a) 


ORIGINAL' PAGE IS 
OF. POOR QUALITY 


log d = 0.9793 log CH -0.1536 (3a) 

Z 0 = 0.164 CH - 0.021 (3b) 

The net radiation absorbed by the canopy R n (W m" 2 ) is calculated according to Rosenthal et. 
al. (1977) using empirical relations for canopy absorption (Kanemasu et. al., 1976) and the net 
radiation above the canopy (Uchijima, 1976) as 

R n = (0.82 S - 90) (1 - exp (-0.398 LAI) ) (4) 

where S is the global insolation (W m~ 2 ) and LAI is the leaf area index. Since in using this empirical 
equation for net radiation I am not calculating the long-wave radiation balance explicitly, the depen- 
dence of net radiation on the canopy temperature is neglected. 

The right hand side of Equation (1) is the resistance analogue for transpiration (van den Honert, 
1948; Gardner, 1964); the numerator is potential difference between leaf (or xylem) and soil water, 
and the denominator is the total resistance for water flow through soil and plant. When the atmo- 
spheric evaporative demand (i.e., the left hand side of Equation (1 )) is zero (which would occur gen- 
erally at night), the leaf water potential i (m) will tend to equilibriate with the soil water potential 
i p s (m) (cf., Biscoe et al., 1976). Boyer (1968), however, notes that because of plant growth the 
leaf water potential before sunrise could be 10 to 20 m lower than the soil water potential. As the 
evaporative demand increases the plant increases the sunction for soil water by decreasing the leaf 
water potential. 

The resistance for water flow through soil R $ is calculated according to Soer (1980) for 
silt loam soil parameters (Clapp and Hornberger, 1978) as 

R s = 550 (<P S / ^t) 2 - 57 (5) 

where the ’saturation* suction is -0.8 m. 
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canopy resistance follows from Equation (7), and the left <$Fth£Q^Rt of Equation 

Ot^ TUer TiuC.H-'T 

(1) gives the rate of transpiration (E). Then from energy balance (cf., Jackson et, al., 1981) 
the canopy-air temperature difference 6T is obtained as 

5T = (R n - L v E) (8) 

C„ p a 

Note that the dependence of net radiation on the canopy temperature is not accounted for explicitly 
due to my use of an empirical equation (4). 

INPUT DATA 

* 

Significant linear relationships between 5T and VPD observed by Idso et. al. (1981) in- 
cluded data from different geographic locations and covered a wide range of weather variables. 
Since weather variables are not totally uncorrelated, the present simulation is done using ob- 
served (U.S. Department of Commerce) hourly weather data for clear skies during summer at 
Phoenix (Az.) and at Baltimore (Md.), A total of sixteen days data consisting of four days for 
Baltimore during June and July and twelve days for Phoenix during May and June. The weath- 
er variables during 1000 to 1400 LST (local standard time) covered a wide range; VPD from 
1.6 to 8.3 kPa. wind speed from 1.5 to 6.5 m/s and air temperature from 23 to 43° C. 

Current weather data do not include global insolation. A model developed by Choudhury 
(1982) was shown to be accurate to 98% or better when compared with careful pyranometer 
observations for clear skies. Considering that the errors in using the empirical equations, such 
as Equation (4), are generally 10 to 15%, the present simulation are done using the insolation values 
calculated from the model. 

The crop height is assumed 0.9 m and the leaf area index 4 as representative for a mature 
soybean crop (Luxmoore et. al.. 1971: Sivakumar et. al.. 1977; Bailey and Davies, 1981a). The 
rooting density is somewhat dependent upon the soil texture and crop variety. In sandy soils 
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the root length per unit area (Lj) is about 9 x IO 3 nT 1 (Robertson et. al., 1980; Jones et. ah, 
1982), and in the surface 0.8 m profile consisting of loam, silt loam and sandy loam (overlying 
several meters of coarse sand and gravel) Arya et. al. (1975) found an L T value of about 1.3 x 
10 4 m" 1 . According to Mayaki et. al. (1976) about 90% of the root dry matter of mature soy- 
beans in silt loam soils is in the surface 0.9 m layer. Raper and Barber (1970) found the root 
length per plant, also in the surface 0.8 m layer, differed among the varieties ranging from 600 
to 800 m. Since the root length per plant of 680 m observed by Arya et. al. (1975) is within 
the range observed by Raper and Barber (1970), an Lj. value of 1.4 x 10 4 m" 1 is assumed to 
be representative for mature soybean crops in silt loam soils. Recognizing the uncertainty in 
and the dynamic nature of root distribution some results will also be given for Lp values of 
(1.4 ± 0.2) x 10 4 m" 1 . The chosen uncertainty in the Lp value reflects the variation in the 
root length observed by Raper and Barber (1970) when compared with Arya et. al. (1975). 

Note that for unirrigated plants Lp would be smaller since Sivakumar et. al. (1977) observed 
that the root length at any depth in silt loam soil decreased when the soil water potential de- 
creased below -20 m. 

SIMULATION RESULTS 

At the soil water potential of -1 m the simulated dependence of 5T on VPD is shown in 
Figure 1 together with the regression line calculated by Idso et. al. (1981) using observations 
on non-stress soybeans at Kansas, Nebraska and North Dakota. In the observations VPD did 
not exceed 6 kPa, and for this range of VPD the simulation generally agrees with the predic- 
tion based on the regression line to within 0,5 K. The data plotted in Idso et. al. indicates 
a standard error of about 1 K with respect to the regression line. Whereas the observed data 
are scattered evenly with respect to the regression line, the simulated results are seen to be 
somewhat biased (lower 5T values), particularly in the VPD range of 2.5 to 4.5 kPa. In fact, a closer 
look at the simulated results indicates a quasi-linear relationship between 5T and VPD, similar to the 


6 


ORIGINAL PAGE IS 
OF POOR QUALITY 


locus for lower bound of 5T calculated by Jackson et al. (1981) for constant stomatal resistance. 
Considering the range of weather variables used in this simulation, the scatter in Figure l supports 
Idso et ai. (1981) in showing that air and dew point temperatures are the significant weather vari- 
ables affecting the canopy temperature during most of clear sky day period. Note that air and dew 
point temperatures also appeared as the significant weather variables affecting evapotranspiration 
from forests (Luxmoore et al., 1981). 

Canopy temperature reflects the integrated response of plant to the existing weather and soil 

water conditions. If the effect of weather variables can be accounted for in terms of air and dew 

point temperatures, e.g. by Figure 1, it may be possible to infer the soil water conditions from the 

canopy temperature measurements, Figure 2 shows the relationship between th'i! 5T values at the 

soil water potential of -1 m (as shown in Figure 1 ) and the 5T values at the potentials of -41, 

-61 and -81 m. Significant relationships between these sets of 5Ts seen in this figure appear poten- 

# 

tially useful in remote sensing of root zone soil water conditions. Referring to Figure 1, let C ' 
assume that the observed value of 5T is 5T 0 . If the soil water potential were -l m (which is close 
to the field capacity) the observed value would have been 5T e (one may calculate this knowing air 
and dew point temperatures using the regression equation of Idso et al.). This 5T e value can now be 
used as the X-coordinate and 5T 0 value as the Y-coordinate in Figure 2 to biacket the soil water 
potential. In the example shown the observed 5T corresponds to the soil water potential of about 
-40 m. 


The highly correlative nature of 5Ts at two different soil water potentials (as seen in Fig. 

2) can be expressed mathematically. If, following Jackson et. al. (1981), we define a crop water 
stress index (CWSI) as 

^ A + 7 (1 + r°/r a ) 


CWSI 


A + 7 (1 + r c .r a ) 


(9) 
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(where r° and r c are, respectively, the canopy resistances at the soil water potentials of -1 m 
and t// s ) then by a simple manipulation of eqns, (l) and (8) one can write 

R n r a 

6>T = ( ) CWSI + (1 - CWSI) 5T° (10) 

C P 

where 6T° and 5T are, respectively, the 5T values at the soil water potentials of -1 m and i// s> 
Changes in soil water potential affects ST through the CWSI. At any soil water potential a 
linear regression analysis between ST and 5T° (for the same VPD) should give a high correlation 
coefficient, and the scatter in this relationship would be due to variable wind conditions. For 
the silt loam soil studied here, a soil water potential of -40 m would correspond to a depletion 
of the plant available soil water by about 75%. From the slope of ST vs. 5T° curve at \p s * 
-41 m shown in Figure 2 we calculate the CWSI - 0.2. Thus, for irrigation scheduling, based 
on 70 to 80% depletion it would be undesirable f ... the CWSI to exceed much beyond 0.2. 

Kramer (1969) suggested that measurements of stomatal resistance would be a direct ap- 
proach to inferring the root zone soil water potential. Brady et. al. (1975) found a second- 
order polynomial relationship between the stomatal resistance (sampling period 1300 to 1500 
CDT) of soybean and soil water potential (silt loam soil), and Sivakumar and Shaw (1978) found 
that the daily average value of the leaf conductance is directly proportional to the soil water 
potential. 

The simulated leaf stomatal resistances from Equation (7) at the soil water potential of -1 m 
are shown in Figure 3 as a function of VPD. In agreement with the observations of Tan and Black 
(1976) the resistance increases as VPD increases. A linear relationship with the VPD would be a 
good approximation except for cases of high evaporative demand (wind speed exceeding 4 m/s 
with VPD greater than about 5 kPa). The soybean canopy temperature data of Kansas shown in 
Idso er. al. (1981) indicate that VPD generally does not exceed 6 kPa. and Figure 3 shows that for 
VPD up to 6 kPa the range of stomatal resistance is 136 to 215 s/m. This range for simulated 
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resistances is in good agreement with the observed range, from 140 to 190 s itt’* 1 , by Brady 
et al. (1975). 

The range of stomatal resistances from 2.2 to 6 kPa at various soil water potentials is shown in 
Figure 4. The increase in stomatal resistance with decreasing soil water potential is consistent with 
observations of Sionit and Kramer (1976). The regression equations calculated by Brady et al. 
(1975) and Sivakumar and Shaw (1978) for soybean crops are also plotted to show agreement with 
the simulation. Note that the range of resistance increases as the soil water potential decreases, as 
was observed by Brady et al. (1975) for soybean and Tan and Black (1976) for forests, Further- 
more, it is seen that the upper bound of the resistance (VPD about 6 kPa) increases more rapidly 
compared to the lower bound (VPD about 2.2 fcJPjj) as the soil water potential decreases. Indeed, 
Tan and Black (1976) observed this dependence of stomatal resistance on the VPD as the soil water 
potential decreased, which led them to conclude that if the measurements of stomatal resistance are 
to be used for inferring soil water conditions, the range of VPD should be constrained to small 
ranges and be prescribed. Clearly, stomatal resistances at high VPDs (say VPD greater than 4 kPa) 
will be more sensitive to the changes in soil water potential, but the simulation shows that such 
resistances are also more sensitive to changes in the wind speed. Therefore, if VPD is to be the only 
weather variable that is prescribed in the stomatal resistance-soil water potential relationship then it 
may be useful to constraint this variable to some moderate value of say 4 kPa. 

The dependence of leaf water potential on VPD during 1300 to 1500 LST is shown in Figure 3 
at the soil water potential of -1 m. and the range of the potential corresponding to stomatal resis- 
tances Figure 4 is shown in Figure 5. The regression equation for soybean calculated by Brady et al. 
(1974) is plotted in Figure 5 to show that the observed data as expressed through the regression 
equation and its standard error or estimate (8.7 m) is fairly well represented in the simulation. Note 
that in contrast to the stomatal resistances shown in Figure 5. the range of leaf water potential 
decreases as the soil water potential decreases. From leaf- and soil-water potential measurements 
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Brady et al, (1974) had concluded that in response to the daily evaporative demand the leaf water 
potential of soybean decreases on the average of 90 to 100 m. In the simulation this decrease is 
essentially given by the difference of leaf and soil water potentials, and Figure 5 shows the magni- 
tude of this decrease to be 70 to 115 m at the soil water potential close to Field capacity, and 60 to 
90 m at the soil water potential of -80 m, In this regard the present simulation agrees marginally 
with Brady et al. (1974). In assessing these comparisons it is pertinent to consider the measurement 
uncertainties. The leaf and soil water potentials measured by Brady et al. (1974) before sunrise did 
not follow a l : 1 relationship which, as such, disagrees with the assumption of plant- and soil-water 
equilibrium. The pressure chamber method used for leaf water potential measurements did not 
include xylem osmotic potential, which Brady et al. estimate could be as much as -10 m when the 
pressure potential is -80 m. If the leaf water potentials calculated from the regression equation of 
Brady et al. (1974) are to be uniformly decreased by 10 m, the prediction based on the regression 
equation would still pass through the simulated range, however, the standard error of estimate for 
the regression equation would go beyond the lower bound of the range for most soil water poten- 
tials. Thus, if the apparent contradiction of data with soil- and plant-water equilibriums is to be 
reconciled by decreasing the measured leaf water potentials by 10 m or more, then the simulation 
will be less satisfactory with respect to these potentials. Carlson et al. (1979) observed that the leaf 
water potentials of two soybean cultivars giving the same stomatal resistance sometimes differed by 
more than 20 m. Observations of Sionit and Kramer (1976) show that the stomatal resistance at 
different growth stages of soybeans sometimes differing by a factor of two or more at the same leaf 
water potential. Therefore, disagreements with respect to the leaf water potential (or the stomatal 
resistance) would indicate the need for a better relationship between the stomatal resistance and 
leaf water potential (Equation 7). The seasonal average predawn and midday values of irrigated soy- 
bean leaf water potential (planted in silt loam soil) observed by Jung and Scott (1980) were, respec- 
tively, -42 and -116 m. If under the assumption of soil- and plant-water equilibrium, the predawn 
value of leaf water potential is assumed to represent the soil water potential then the observed 
midday value of the leaf water potential (-116 m) would not be inconsistent with the present simu- 
lation. However, the corresponding observed seasonal average midday stomatal resistance of 60 s/m 
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is about a factor of three smaller than that predicted by the simulation. Thus, although the present 
simulation provides some qualitative understanding of a few independent sets of observations, it also 
shows the need for a better quantification of the stomatal resistance - leaf water potential relation- 
ship for different soybean cultivars. It is also pertinent to quantify the mechanics of osmotic regula- 
tion and the role of leaf turgor in stomatal control (cr., Turner et al., 1978; Zur and Jones, 1981). 

The above simulation results are obtained for root length per unit area (L T ) value of 1.4 x 
10 4 m -1 . This plant characteristic vary with time, and Figure 6 illustrates the 8T vs. VPD relation- 
ship for L t values of 1.2 x 10 4 and 1.6 x 10 4 m“ l . It is seen that the 5T values decrease as the L T 
value increases, i.e., as the rooting density of a plant increases the canopy temperature decreases. 

The dependence of 5T on L T ere not the same at all VPDs; for example, at the VPD of 2 kPa the 
change in 5Ts is about 0.3 K, and at the VPD of 7 kPa the change in 5Ts is about 1.0 K. The simula- 
tion suggests that some scatter in the observations of Idso et al. (1981) may also have arised from 
differing rqotins;* tfsis&kies. Also, since the observed rooting density in sandy soils (Robertson et al., 
1980; Jones et al., 1981) appears to be lower than in loam soils (Arya et al., 1975; Raperand Barber, 
1970), the simulation would suggest that the canopy temperature of soybeans planted in sandy soil 
to be somewhat higher compared to the crop planted in loam soils. 

SUMMARY 

A soil-plant-atmosphere model for soybean was described and simulation results using observed 
hourly weather data were compared with observed canopy temperatures, stomatal resistances, and 
leaf water potentials, The air and dew point temperatures are found to be the significant weather 
variables affecting the canopy temperatures. Under identical weather conditions, the model gives a 
lower canopy temperature for a soybean crop with a higher rooting density. A knowledge of crop 
rooting density, in addition to air and dew point temperatures, is needed in interpreting infrared 
radiometric observations for soil water status. The observed dependence of stomatal resistance on 
VPD and soil water potential is fairly well reproduced. Analysis of the simulated leaf water poten- 
tials indicated overestimation, possibly due to differences in the cultivars. 
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Figure 1. Simulated variation of soybean canopy-air temperature difference with air 
vapor pressure deficit. The regression equation calculated by Idso et. al. 
(1981) from observations on unstress canopies is plotted. The soil water 
potential is - 1 m. 


17 






LEAF STOMATAL RESISTANCE {s/m) LEAF WATER POTENTIAL (m) 


OmGINAL P flGE ig 
° F P00R QUALITY 



60 




RANGE OF LEAF STOMATAL RESISTANCES (s/m> 




SOI L WATER POTENTIAL (m) 

Figure 4. The simulated range of leaf stomatal resistances at various soil water po- 
tentials. The dependencies inferred from observations bv Brady et. al. 
(1975) and Sivakumar and Shaw (1978) are also plotted. 
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Figure 5. The range of leaf water potentials at various soil water potentials. The re- 
gression equation calculated by Brady et. ai. (1974) is shown. 
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Figure 6. Hie effect of differing rooting densities on the canopy-air temperature difference 
at the soil water potential of -1 m. The canopy temperature decreases as the root 
length per unit area (L T ) increases. The rooting density affects the canopy temper- 
atures more at higher vapor pressure deficits. 
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Simulated variation of soybean canopy-air temperature difference with air vapor 
pressure deficit. The regression equation calculated by Idso et. al. (1981) from 
observations on unstress canopies is plotted. The soil water potential is - 1 m. 

Relationship of the canopy-air temperature differences at the soil water potentials 
of - 41, - 61, and -81m with respect to the soil water potential of - 1 m. With 
decreasing soil water potential the scatter increases due to differing wind speeds. 

Simulated dependence of leaf stomatal resistance and leaf water potential on the 
air vapor pressure deficit. The soil water potential is •* l m. Straight lines are 
drawn to show the trend. 

The simulated range of leaf stomatal resistances at various soil water potentials. 

The dependencies inferred from observations by Brady et. al. (1975) and Sivakumar 
and Shaw (1978) are also plotted. 

The range of leaf water potentials at various soil water potentials. The regression 
equation calculated by Brady et. al. (1974) is shown. 

The effect of differing rooting densities on the canopy-air temperature difference 
at the soil water potential of - 1 m. The canopy temperature decreases as the root 
length per unit area (L T ) increases. The rooting density affects the canopy temper- 
atures more at higher vapor pressure deficits. 


